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We place constraints on the number of relativistic degrees of freedom and on the baryon asymmetry 
at the epoch of Big Bang Nucleosynthesis (BBN) and at recombination, using cosmic background 
radiation (CBR) data from the Wilkinson Microwave Anisotropy Probe (WMAP), complemented 
by the Hubble Space Telescope (HST) Key Project measurement of the Hubble constant, along 
with the latest compilation of deuterium abundances and H n region measurements of the primordial 
helium abundance. The agreement between the derived values of these key cosmological and particle 
physics parameters at these widely separated (in time or redshift) epochs is remarkable. From the 
combination of CBR and BBN data, we find the 2a ranges for the effective number of neutrinos 
N v and for the baryon asymmetry (baryon to photon number ratio n) to be 1.7-3.0 and 5.53-6.76 
xlO -10 , respectively. 
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I. INTRODUCTION 

The concordance model of cosmology, with dark en- 
ergy, dark matter, baryons, and three flavors of light neu- 
trinos, provides a consistent description of BBN (~ 20 
minutes), the CBR (~ 380 Kyr), and the galaxy forma- 
tion epochs of the universe (> 1 Gyr). The standard 
model has received recent confirmation from the WMAP 
precision measurements of the CBR temperature and po- 
larization anisotropy spectra [1]. However, despite the 
impressive successes of the standard model in describing 
a wide range of cosmological data, the possibility remains 
that there could be non-standard model contributions to 
the total energy density in the radiation era from addi- 
tional relativistic particles. 

In this paper new constraints are placed on any physics 
beyond the standard model that contributes to the en- 
ergy density like radiation (i.e., decreases with the expan- 
sion of the universe as the fourth power of the scale factor, 
independent of the sign of that contribution) . While such 
new physics may or may not be due to extra relativistic 
degrees of freedom, it is assumed that the non-standard 
contribution to the energy density may be parameterized 
as such. Simultaneously, constraints are placed on the 
baryon density at widely differing epochs in the evolution 
of the universe. The keys to these constraints are the 
recently released measurements of the CBR anisotropy 
spectra by the WMAP collaboration, the most recent 
compilation of high redshift, low metallicity deuterium 
abundances [2] and 4 He abundances relevant to BBN. 



energy density receives its dominant contribution from 
all the relativistic particles present (the evolution of the 
universe is said to be "radiation-dominated" (RD)). In 
the standard cosmology, prior to annihilation, these 
relativistic particles are: photons, e pairs and three fla- 
vors of left-handed (i.e., one helicity state) neutrinos (and 
their right-handed antineutrinos). Then, the energy den- 
sity is 



PTOT = PR. = P-f + Pe + 'iPu 
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where p 1 is the energy density in photons (which by to- 
day have redshifted to become the CBR photons at a 
temperature of about 2.7 K). 

In "standard" BBN (SBBN) it is assumed that the neu- 
trinos are fully decoupled prior to annihilation and do 
not share in the energy transferred from the annihilating 
e ± pairs to the CBR photons. In this approximation, the 
photons in the post-e^ annihilation universe are hotter 
than the neutrinos by a factor T 7 /T„ = (11/4) 1 / 3 , and 
the relativistic energy density is 



Pn, = p-y + ?>Pv = 1.6813^ 



(2) 



During the RD epoch the age and the energy density 
are related by |/9R,i 2 = 1 (we have chosen units in which 
8nG = 1), so that once the particle content (pn) is speci- 
fied, the age of the universe is known as a function of the 
CBR temperature. In the standard model, 

Pre - e ± annihilation : t T 2 = 0.738 MeV 2 s, (3) 



II. MODIFIED RELATIVISTIC ENERGY 
DENSITY 



Post 



2^ annihilation 



t T 2 = 1.32 MeV 2 s. 



(4) 



The cosmology of interest here begins when the uni- 
verse is already a few tenths of a second old and the tem- 
perature is a few MeV. At such early epochs the total 



The most straightforward variation of the standard 
cosmology is "extra" energy contributed by new, light 
(relativistic at BBN) particles "V. These might, but 
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need not be sterile neutrinos. When the X are decou- 
pled, in the sense that they don't share in the energy 
released by annihilation, it is convenient to account 
for the extra contribution to the standard-model energy 
density by normalizing it to that of an "equivalent" neu- 
trino [3], 



III. 



CONSTRAINTS ON N v FROM THE CBR 



Px = AN vPu = -AN vPl . 



(5) 



For SBBN AN V = 0, where AiV,, = 3 + N v . For each ad- 
ditional "neutrino-like" particle (i.e., any two-component 
fermion), if Tx = T v , then AN V — 1; if X is a scalar, 
AA„ = 4/7. However, it may well be that the X have 
decoupled even earlier in the evolution of the universe 
and have failed to profit from the heating when vari- 
ous other particle-antiparticle pairs annihilated (or un- 
stable particles decayed). In this case, the contribution 
to AN U from each such particle will be < 1 (< 4/7). 
We emphasize that, in principle, we are considering any 
term in the energy density which scales like a -4 , where 
a is the scale factor. In this sense, the modification to 
the usual Friedman equation due to higher dimensional 
effects, as in the Randall-Sundrum model [4] (see also, 
[5, 6, 7, 8, 9, 10, 11, 12]), can be included as well. An 
important interest in this latter case is that it permits 
the possibility of a negative contribution to the radiation 
density (AN U < 0; N v < 3). 

In the presence of such a modification to the relativistic 
energy density, the pre-e* annihilation energy density in 
Eq. (1) is changed to, 



(PR.) 



43 



7AN„ 
43 



Pi- 



(6) 



Any extra energy density (AiV„ > 0) speeds up the ex- 
pansion of the universe so that the right-hand side of the 
time-temperature relation in Eq. (3) is smaller by the 
square root of the factor in parentheses in Eq. (6), 



(t/t') p 

r 



7AJV, 
43 



1/2 



(1 + 0.1628AiV„) */ 2 , (7) 



where t' is the age of the universe with the extra energy 
density. In the post-e 1 * 1 annihilation universe the extra 
energy density is diluted by the heating of the photons, 
so that, 



(PR 



post 



1.6813 (l + 0.1351A^)p 7 , 



(8) 



and 



S post = (t/t') post = (1 + 0.1351A7V,) 1 / 2 . (9) 



This latter expression (Eq. 9) is also relevant for the mod- 
ification to the spectrum of temperature fluctuations in 
the CBR (when compared with the standard N v = 3 
case). 



The competition between gravitational potential and 
pressure gradients is responsible for the peaks and 
troughs in the CBR power spectrum. The redshift of 
matter-radiation equality, 



z e9 = 2.4xl0 4 ^ 

^post 



(10) 



affects the time (redshift) duration over which this com- 
petition occurs. Here, lum = Vtuih 2 is the total mat- 
ter density (comprised, for nearly massless neutrinos, 
of baryons and cold dark matter) and h (Ho = 100/z 
km/s/Mpc) is the normalized Hubble constant. The di- 
rect correlation between u>m and AiV„ is evident in Fig. 1 
which results from our analysis described below. The 
primary effects of relativistic degrees of freedom (other 
than photons) on the CBR power spectrum result es- 
sentially from changing the redshift of matter-radiation 
equality. If the radiation content is increased, matter- 
radiation equality is delayed, and occurs closer (in time 
and/or redshift) to the epoch of recombination. 



< 




0.4 



FIG. 1: The CBR degeneracy between u>m and AiV„ is evident 
from the la and 2a contours from the WMAP data. 



The redshift of matter-radiation equality is important 
for two reasons [13]: 

• Radiation causes potential decay which blueshifts 
the photons because they do not have to climb out 
of such deep wells. Moreover, the concurrent de- 
cay in the spatial curvature doubles the blueshift 
effect by contracting the wavelength of the photons 
relative to the pure cosmological expansion. 

• In the matter dominated (MD) era before recombi- 
nation, the density contrast (Sp/p) of the pressure- 
less cold dark matter (CDM) grows unimpeded (as 
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i 2 / 3 ) while the density contrast of the baryons is 
either oscillating or decaying. The longer this pre- 
recombination MD era lasts, the more suppressed 
are the amplitudes of the peaks. 

Conversely, if matter-radiation equality is delayed, the 
gravitational potential is dominated by the photon- 
baryon fluid closer to recombination resulting in a more 
pronounced peak structure. 

An increase in the relativistic content causes the uni- 
verse to be younger at recombination with a correspond- 
ingly smaller sound horizon s*. Since the location of the 
n th peak scales roughly as nnD^/s^ (where £>* is the co- 
moving angular diameter distance to recombination), the 
peaks shift to smaller angular scales (larger I) and with 
greater separation. These features are clearly visible in 
Fig. 2. 

The heights and locations of the peaks also depend on 
the history of the universe after recombination. At the 
end of matter domination and the onset of dark energy 
domination, further and much slower (compared to that 
in the radiation epoch) potential decay occurs. The more 
gradual potential decay causes the induced anisotropy to 
be suppressed by a factor of I. The amplification of the 
power in the lowest I's from this late decay serves as a 
probe of dark energy (or another probe of the matter 
content in a flat universe). In principle, the degeneracy 
between AN U and ujm is broken by this effect and by 
the accompanying change in the redshift at which the 
matter dominated epoch ends. However, note that the 
lowest multipoles also have the largest cosmic variance. 

7000 I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — i — I 




Multipole 1 

FIG. 2: The power spectrum for the best-fit (N v = 2.75) to 
the WMAP data is the solid line. With all other parame- 
ters and the overall normalization of the primordial spectrum 
fixed, the spectra for N v — 1, N v — 5 and N v = 7 are the 
dotted, dot-dashed and dashed lines, respectively. The data 
points represent the binned TT power spectrum from WMAP. 



The TT and TE power spectra are computed using 



the Code for Anisotropies in the Microwave Background 
or CAMB [14] which is a parallelized version of CMB- 
FAST [15]. The Universe is assumed to be flat, in accord 
with the predictions of inflation [16], and the dark en- 
ergy is assumed to behave as a cosmological constant A. 
The restriction of a flat geometry allows us to relate the 
dark energy and matter densities at the present time: 
SIa = 1 — Qai- The angular power spectrum is calculated 
on a grid defined by h, the baryon density wb = Sl-ah 2 (or 
rj W = 10 10 rtB/^ 7 = 274<x>b)j wm, the number of equiva- 
lent neutrinos N u , the reionization optical depth r, and 
the spectral index n s of the primordial power spectrum. 
Two priors are imposed to largely break the degener- 
acy between lom and AN V . For h a top-hat distribu- 
tion is chosen corresponding to the HST measurement, 
h = 0.72 ± 0.08 [17], and we require that the universe be 
older than the globular clusters (which, at 2<r, are older 
than 11 Gyr [18]). For comparison, we also consider the 
case when the age of the universe t exceeds 12 Gyr. 
Our top-hat grid, consisting of over 10 million points, 

is: 

• 0.64 < h < 0.8 in steps of size 0.02. 

• 0.018 <w B < 0.028 in steps of size 0.001. 

• 0.11 < ojm < 0.27 in steps of size 0.01 and ujm = 
0.07, 0.3. 

• 1 < N v < 3.5 in steps of size 0.25, 4 < N v < 9 in 
steps of size 0.5 and N v = 0, 0.5. 

• < r < 0.3 in steps of size 0.025. 

• 0.90 <n s < 1.02 in steps of size 0.01 and n s = 0.80, 
0.84, 0.88, 1.04, 1.08, 1.12, 1.16, 1.20. 

• The normalization of the spectrum is a continuous 
parameter. 

The first year WMAP data are in the form of 899 
measurements of the TT power spectrum from I = 2 
to I = 900 [19] and 449 data points of the TE power 
spectrum [20] . Although the effect of relativistic degrees 
of freedom on the TE spectrum is insignificant, it is in- 
cluded in our analysis for completeness. The likelihood of 
each model of our grid is computed using Version 1.1 of 
the code provided by the WMAP collaboration [21]. The 
code computes the covariance matrix under the assump- 
tion that the off-diagonal terms are subdominant. This 
approximation breaks down for unrealistically small am- 
plitudes. When the height of the first peak is below 5000 
/iK 2 (which is many standard deviations away from the 
data), only the diagonal terms of the covariance matrix 
are used to compute the likelihood. 

The best-fit parameters are h = 0.68, lub = 0.023 
(?7io = 6.3), lom = 0.14, N v = 2.75, r = 0.13, and 
n s = 0.97 with a x 2 = 1429.13 for 1341 degrees of free- 
dom. The allowed parameter space in the Tyio-AA^ plane 
is shown in Fig. 3. The solid (dotted) lines correspond to 
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FIG. 3: The la and 2a contours in the ^lo-AiV,, plane from 
WMAP data. The solid (dotted) lines correspond to to > 11 
(12) Gyr. The cross marks the best-fit at lub = 0.023 and 
AN V = -0.25. 



the la and 2a regions 1 for to > H (12) Gyr. The cross 
identifies the best-fit point. Note that while this best fit 
point lies at N v < 3, the N„ distribution is very broad. 
After marginalizing over 77, the 2a range in N v extends 
from 0.9 {AN V = -2.1) to 8.3 (AN„ = 5.3). Although 
similar CBR analyses (see Ref. [22]) have included differ- 
ent, additional data to that from WMAP alone, making 
direct comparisons difficult, our results are in good agree- 
ment with them. The prior on to has a significant effect 
on the allowed values of AN V [23] for a simple reason. 
Since flatness is assumed, to depends only on the matter 
content and the Hubble parameter via 



Ho to — - 



1 



3 y/T^Ti 



:hl 



M 



1 + y/l - O.M 



(11) 



The combination of the HST prior on h and the to prior 
restricts lom and help to break the degeneracy between 
AN V and wm- 

The best fit WMAP-determined baryon density is 
7710 = 6.30 (wb = 0.0230), in excellent agreement with 
Spergel et at [24] and other similar analyses [22]. The 
CBR 2a range extends from 7710 = 5.58 (lob = 0.0204) to 
7710 = 7.26 (w B = 0.0265). 

These CBR constraints on N v and lob apply to epochs 
in the evolution of the universe 380 Kyr. An impor- 
tant test of the standard models of cosmology and parti- 



1 For 2-dimensional constraints, the lcr, 2a and 3(7 regions are 
defined by A\ 2 = 2.3, 6.17 and 11.83, respectively. 



cle physics is to compare them with corresponding con- 
straints from the much earlier epoch probed by BBN. 



IV. THE ROLES OF N v AND 7710 IN BBN 

At T ~ few MeV, the neutrinos are beginning to de- 
couple from the 7 — plasma and the neutron to proton 
ratio, crucial for the production of primordial 4 He, is de- 
creasing. As the temperature drops below ~ 2 MeV, 
the two-body collisions between neutrinos and e ± pairs, 
responsible for keeping the neutrinos in thermal equilib- 
rium with the electron-positron - photon plasma become 
slow compared to the universal expansion rate and the 
neutrinos decouple, although they do continue to interact 
with the neutrons and protons via the charged-current 
weak interactions. Prior to e ± annihilation, when the 
temperature drops below ~ 0.8 MeV and the universe is 
~ 1 second old, these interactions, interconverting neu- 
trons and protons, become too slow (compared to the uni- 
versal expansion rate) to maintain n — p equilibrium and 
the neutron-to-proton ratio begins to deviate from (ex- 
ceeds) its equilibrium value ((n/p) eq = exp(— Am/T)), 
where Am is the neutron-proton mass difference. Beyond 
this point, often described as neutron-proton "freeze- 
out", the n/p ratio continues to decrease, albeit more 
slowly than would have been the case in equilibrium. 
Since there are several billion CBR photons for every 
nucleon (baryon) , the abundances of any complex nuclei 
are entirely negligible at these early times. 

We note here that if there is an asymmetry between 
the numbers of v e and v e ("neutrino degeneracy"), de- 
scribed by a chemical potential fi e , then the equilib- 
rium neutron-to-proton ratio is modified to (n/p) eq — 
exp(— Am/T — /i e /T). In place of the neutrino chemical 
potential, it is convenient to introduce the dimensionlcss 
degeneracy parameter £ e = fi e /T. A positive chemical 
potential (£ e > 0; more v e than v e ) leads to fewer neu- 
trons and less 4 He will be synthesized in BBN. 

BBN begins in earnest after e annihilation, at T w 
0.08 MeV (t«3 minutes), when the number density of 
those CBR photons with sufficient energy to photodis- 
sociate deuterium (those in the tail of the black body 
distribution) is comparable to the baryon density. By 
this time the n/p ratio has further decreased (the two- 
body reactions interconverting neutrons and protons hav- 
ing been somewhat augmented by ordinary beta decay; 
t„ = 885.7 sec), limiting (mainly) the amount of helium- 
4 which can be synthesized. As a result, the predictions 
of the primordial abundance of 4 He depend sensitively 
on the early expansion rate and on the amount - if any 
- of a v e — D e asymmetry. 

In contrast to 4 He, the BBN-predicted abundances of 
deuterium, helium-3 and lithium- 7 (the most abundant of 
the nuclides synthesized during BBN) are determined by 
the competition between the various two-body produc- 
tion/destruction rates and the universal expansion rate. 
As a result, the D, 3 He, and 7 Li abundances are sensitive 
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to the post-e annihilation expansion rate, while that of 
4 He depends on both the pre- and post-e 1 * 1 annihilation 
expansion rates; the former determines the "freeze-in" 
and the latter modulates the importance of beta decay 
(see, e.g., Kneller & Steigman [25]). Also, the primordial 
abundances of D, 3 He, and Li, while not entirely insensi- 
tive to neutrino degeneracy, are much less effected by a 
non-zero £ e (e.g., [26]). 

Of course, the BBN abundances do depend on the 
baryon density which fixes the nuclear reactions rates and 
also, through the ratio of baryons to photons, regulates 
the time/temperature at which BBN begins. As a result, 
the abundances of at least two different relic nuclei are 
needed to break the degeneracy between the baryon den- 
sity and a possible non-standard expansion rate resulting 
from new physics or cosmology, and/or a neutrino asym- 
metry. In this paper only the former possibility is consid- 
ered; in another publication several of us (along with P. 
Langacker) have explored the consequences of neutrino 
degeneracy and we studied the modifications to the con- 
straints on AN V when both of these non-standard effects 
are simultaneously included. 



weakly, correlated with AN V . In contrast to D (and to 
3 He and Li), since the incorporation of all available neu- 
trons into He is not limited by the very rapid nuclear 
reaction rates, the 4 He mass fraction is relatively insen- 
sitive to the baryon density, but it is very sensitive to 
both the pre- and post-e 1 * 1 annihilation expansion rates 
(which control the neutron-to-proton ratio). The faster 
the universe expands, the more neutrons are available 
for 4 He. The very slow increase of Y with 7710 is a re- 
flection of the fact that for higher baryon density, BBN 
begins earlier, when there are more neutrons. As a re- 
sult of these complementary correlations, the pair of pri- 
mordial abundances yu = 10 5 (D/H) and the 4 He mass 
fraction Y, provide observational constraints on both the 
baryon density and the universal expansion rate when 
the universe was some 20 minutes old. Comparing these 
to constraints when the universe was some 380 Kyr old, 
from the WMAP observations of the CBR spectra, pro- 
vides a test of the consistency of the standard models of 
cosmology and of particle physics and further constrains 
the allowed range of the present baryon density of the 
universe. 
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FIG. 4: Isoabundance curves for D and 4 He in the 7710 - 
AN„ plane. The solid curves are for 4 He (from top to bottom: 
Y = 0.25, 0.24, 0.23). The dashed curves are for D (from left 
to right: 10 5 (D/H) = 3.0, 2.5, 2.0. The data point with error 
bars corresponds to yo = 2.6±0.4 and Y = 0.238±0.005; see 
the text for discussion of these abundance values. 



While the abundances of D, 3 He, and Li are most sen- 
sitive to the baryon density (77), the 4 He mass fraction 
(Y) provides the best probe of the expansion rate. This 
is illustrated in Fig. 4 where, in the AN V - 7710 plane, 
are shown isoabundance contours for D/H and Y (the 
isoabundance curves for 3 He/H and for Li/H, omitted 
for clarity, are similar in behavior to that of D /H) . The 
trends illustrated in Fig. 4 are easy to understand in the 
context of the discussion above. The higher the baryon 
density (7710), the faster primordial D is destroyed, so the 
relic abundance of D is anticorrelated with 7710- But, the 
faster the universe expands (AN V > 0), the less time is 
available for D-destruction, so D/H is positively, albeit 



V. PRIMORDIAL ABUNDANCES 

It is clear from Fig. 4 that while D (and/or 3 He and/or 
7 Li) largely constrains the baryon density and 4 He plays 
the same role for AN V , there is an interplay between 7710 
and AN„ which is quite sensitive to the adopted abun- 
dances. For example, a lower primordial D/H increases 
the BBN- inferred value of 7710, leading to a higher pre- 
dicted primordial 4 He mass fraction. If the primordial 
4 He mass fraction derived from the data is "low", then 
a low upper bound on AA„ will be inferred. There- 
fore, it is crucial to make every effort to avoid biasing 
any conclusions by underestimating the present uncer- 
tainties in the primordial abundances derived from the 
observational data. For this reason deuterium is adopted 
as the baryometer of choice. Primarily, this is because 
its observed abundance should have only decreased since 
BBN [27], but also because the deuterium observed in 
the high redshift, low metallicity QSO absorption line 
systems (QSOALS) should be very nearly primordial. 
In contrast, the post-BBN evolution of 3 He and of 7 Li 
are considerably more complicated, involving competi- 
tion between production, destruction, and survival. As a 
result, at least so far, the current, locally observed (in the 
Galaxy) abundances of these nuclides have been of less 
value in constraining the baryon density than has deu- 
terium. Nonetheless, inferring the primordial D abun- 
dance from the QSOALS has not been without its difficul- 
ties, with some abundance claims having been withdrawn 
or revised. Presently there are 5-6 QSOALS with rea- 
sonably firm deuterium detections [2, 28, 29, 30, 31, 32]. 
However, there is significant dispersion among the abun- 
dances and the data fail to reveal the anticipated "deu- 
terium plateau" at low metallicity or at high redshift [33] . 
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Furthermore, subsequent observations of the D'Odorico 
et al. [32] QSOALS by Levshakov et al. [34] revealed a 
more complex velocity structure and led to a revised - 
and more uncertain - deuterium abundance. This sen- 
sitivity to the often poorly constrained velocity struc- 
ture in the absorbers is also exposed by the analyses of 
published QSOALS data by Levshakov and collabora- 
tors [35, 36, 37], which lead to consistent, but somewhat 
higher deuterium abundances than those inferred from 
"standard" data reduction analyses. In the absence of a 
better motivated choice, here we adopt the five abun- 
dance determinations collected in the recent paper of 
Kirkman et al. [2]. The weighted mean value of j/tj is 
2.6 2 . But, the dispersion among these five data points 
is very large. For this data set x 2 = 15.3 for four degrees 
of freedom, suggesting that one or more of these abun- 
dance determinations may be in error, perhaps affected 
by unidentified and unaccounted for systematic errors. 
For this reason, we follow the approach advocated by 
[30] and [2] and adopt for the uncertainty in yu the dis- 
persion divided by the square root of the number of data 
points. Thus, the primordial abundance of deuterium to 
be used here is chosen to be: yr> = 2.6 ± 0.4. For SBBN 
(iV„ = 3, £ e = 0), at ±lcr this corresponds to a baryon 
density 7710 = 6.1 j^g (lu b = 0.022 ± 0.002) 3 . 

A similar, less than clear situation exists for determi- 
nations of the primordial abundance of 4 He. At present 
there are two, largely independent, estimates based on 
analyses of large data sets of low-metallicity, extragalac- 
tic Hn regions. The "IT" [38, 39] estimate of Y(IT) 
= 0.244± 0.002, and the "OS" determination [40, 41, 42] 
of Y(OS) = 0.234 ± 0.003 which differ by nearly 3a. 
The recent analysis of high quality observations of a rel- 
atively metal-rich (hence, chemically evolved and post- 
primordial) Hn region in the Small Magellanic Cloud 
(SMC) by Peimbert, Peimbert, and Ruiz (PPR) [43] 
yields an abundance Y SM c = 0.2405 ± 0.0018. When 
PPR extrapolated this abundance to zero metallicity, 
they found Y(PPR) = 0.2345 ± 0.0026, lending support 
to the OS value. These comparisons of different observa- 
tions and analyses suggest that unaccounted systematic 
errors may dominate the statistical uncertainties. In- 
deed, Gruenwald, Steigman, and Viegas [44] argue that 
if unseen neutral hydrogen in the ionized helium region 
of the observed H 11 regions is accounted for, the IT es- 
timate of the primordial abundance should be reduced 
to Y(GSV) = 0.238 ± 0.003 (see also [45, 46]). Here, 
we adopt this latter estimate for the central value but, 
as we did with deuterium, the uncertainty is increased 
in an attempt to account for likely systematic errors: Y 



2 This differs from the result quoted in Kirkman et al. because 
they have taken the mean of log(j/D) and then used it to infer 
y D (y D = I0< l °s(vu)>). 

3 We have purposely avoided quoting the baryon density to more 
significant figures than is justified by the accuracy of the D- 
abundancc determination. 



= 0.238±0.005, leading to a 2a range, 0.228 < Y < 0.248; 
this range is in accord with the estimate adopted by 
Olive, Steigman, and Walker (OSW) [47] in their review 
of SBBN. Although we will comment on the modification 
to any conclusions if Y(IT) is substituted for Y(OSW), 
Figs. 4 - 8 are shown for y D = 2.6 ± 0.4 and Y(OSW) = 
0.238 ±0.005. 



VI. STANDARD BBN 
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FIG. 5: The band is the SBBN predicted relation between 
the primordial abundances of D and 4 He, including the errors 
(±lc) in those predictions from the uncertainties in the nu- 
clear and weak interaction rates. The point with error bars is 
for the relic abundances of D and 4 He adopted here (see the 
text). 



Before proceeding to our main goal of constraining new 
physics using BBN, it is worthwhile to set the scene by 
considering the standard model case (N v = 3, £ e = 0) 
first. The result of this comparison is well known: there 
is a "tension" between the primordial abundances of D 
and 4 He inferred from the observational data and those 
predicted by SBBN. For example, if y^ is used to fix the 
baryon density (y D = 2.6 ± 0.4) then, at ±la, r]f^ BN = 
6.lto 7 5 (w B = 0.022±0.002),the corresponding predicted 
4 He abundance is Y = 0.248 ± 0.001 (la), which is some 
2a higher than either the OSW or the IT estimates. This 
is illustrated in Fig. 5 which shows the SBBN-predicted 
relation between the relic abundances of D and 4 He along 
with the Y(OSW) abundance estimate adopted here. 
The D-inferred baryon density is in excellent agreement 
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with the baryon density determined independently (non- 
BBN; N v = 3) by Spergel et al. [24] from a combination 
of CBR and Large Scale Structure data (2dF + Lyman 
a): ry^-BBN = 6 . 14 ± .25 (w B = 0.0224 ± 0.0009) 4 . 
Thus, it appears that 4 He is the problem: the primordial 
abundance of 4 He is smaller than predicted for SBBN 
given either the observed deuterium abundance or the 
non-BBN inferred baryon density. At the same time we 
strongly emphasize that this "discrepancy" is only at the 
~ 2cr level and we should celebrate the excellent agree- 
ment between the baryon density determined when the 
universe was only 20 minutes old and when the universe 
was some 380 Kyr old (CBR). 



VII. NON-STANDARD BBN: AN V 

As noted above, for SBBN (N u = 3) the observation- 
ally inferred primordial abundance of 4 He is too small 
(by ~ 2a) for the baryon density inferred either from the 
D abundance or from the non-BBN analysis of Spergel et 
al. [24]. This suggests that the early universe expansion 
rate may have been too fast, leaving too many neutrons 
available for the synthesis of 4 He. If this tension between 
4 He and D (or, between 4 He and wb) should persist, it 
could be a signal of non-standard physics corresponding 
to S < 1 {AN V < 0). Indeed, in Fig. 4 it can be seen 
that for the adopted primordial abundances of D and 
4 He, there is a "perfect" fit ( X 2 = 0) for AN V w -0.7 
(N u ps 2.3) and 7710 « 5.7. Although N u = 3 is only 
disfavored by ~ 2a, any increase in the early universe 
expansion rate (S > 1, N v > 3) is strongly disfavored. 
This is illustrated in Fig. 6 which shows the la, 2a and 
3a contours in the 7710-AiV^ plane for the adopted D and 
4 He (OSW) abundances. The shape of these contours 
reflects our much discussed complementarity between D 
and He: D provides the best constraint on the baryon 
density while He is most sensitive to the early universe 
expansion rate, the latter providing an excellent probe of 
possible new physics. 

With reference to that Fig. 6 we note that even one 
extra, fully thermalized neutrino (AN V = 1) is strongly 
disfavored. This seemingly eliminates the sterile neutrino 
suggested by the LSND experiment [48] . For the LSND 
parameters, in the absence of significant neutrino asym- 
metry, this "sterile" neutrino would have been mixed 
with the active neutrinos and thermalized prior to neu- 
trino decoupling (prior to BBN) [49]. If such a neutrino 
were to exist, the "new" standard model would corre- 
spond to N v = 4. This would be a disaster since for 
N v = 4 and the OSW 4 He abundance the minimum x 2 



4 It should be noted that from the CBR alone Spergel et al. find 
jjgBR _ g g ± 3 _ o 024 ± 0.001). It is this value which 
is most directly comparable to our CBR and joint BBN/CBR 
results. 
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FIG. 6: The la, 2<r and 3<r contours in the r7io-AY„ plane for 
the adopted D and 4 He (OSW) abundances (solid lines). The 
cross marks the best fit BBN point. The la and 2a contours 
from WMAP (dashed lines) are shown for comparison. 



(which now occurs at 7710 ps 6.1) is greater than 20. The 
situation is even worse for the IT 4 He abundance since 
the smaller uncertainty in Y forces a much smaller baryon 
density (7710 » 4.7), with xLin > 60! 



A. Requiring N v > 3 

Since, as is well known from LEP [50], there are 
three flavors of active, left-handed neutrinos (and their 
right-handed antiparticles) , any extra contributions to 
the relativistic energy density at BBN should result in 
N v > 3. Of course, "new physics" in the form of 
non-minimally coupled fields (e.g., [25, 51] and refer- 
ences therein) or from higher-dimensional phenomena 
{e.g., [4, 5, 6, 7, 8, 9, 10, 11, 12]) may result in an ef- 
fective N u < 3. If, however, the class of non-standard 
physics of interest is restricted to AN V > 0, then the 
BBN constraints presented above (and those from the 
CBR) will change. With a prior of > 3, the best 
fit BBN-determined values of the baryon-to-photon ra- 
tio and AN„ (for Y(OSW)) shift from 7710 ~ 5.7 and 
AN„ « -0.7, to 7710 » 5.9 and AiV„ = 0. The value 
of Xmin changes to 4.2. The corresponding confidence 
contours in the rjiQ-AN^ plane are shown in Fig. 7. 



VIII. JOINT CONSTRAINTS AND SUMMARY 

As may be seen from Figs. 6 and 7, the agreement be- 
tween the values obtained for AN„ and 7710 from WMAP 
and from BBN separately is excellent. Guided by this, 
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FIG. 7: The la, 2a and 3<r contours (solid lines) in the Jjio- 
AiV„ plane for N„ > 3 and the adopted D and 4 He (OSW) 
abundances. The corresponding la and 2a contours from 
WMAP (dashed lines) are shown for comparison. 
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FIG. 8: The la, 2a and 3a contours in the TJio-AiV,, plane 
from a combination of WMAP data and the adopted D and 
4 He (OSW) abundances. 



the BBN and CBR results are combined to obtain the 
joint fit in the r\ - AN V plane shown in Fig. 8. To a 
good first approximation, BBN (and primordial He) de- 
termines AN V while WMAP fixes 7710 (with some help 
from BBN and primordial D). The corresponding figure 
for N v > 3 for the joint fit has not been shown because 
again the AiV„ range is almost identical to that from 
BBN alone (Fig. 7). 

Our results are summarized in Tables I and II and 





N v (2a range) 


7/10 (2a range) 


T 1 T" TV T A TTl 

WMAP 


0.9 - 8.3 


5.58 - 7.26 


yo + Y(OSW) 


1.7 - 3.0 


4.84 - 7.11 


Jto + Y(IT) 


2.4 - 3.0 


5.06 - 7.33 


WMAP + y D + Y(OSW) 


1.7 - 3.0 


5.53 - 6.76 


WMAP + y D + Y(IT) 


2.4 - 3.0 


5.58 - 6.71 



TABLE I: The 2a ranges (for 1 degree of freedom) of N v and 
7710 from analyses of WMAP data, deuterium and helium 
abundances and their combinations. The WMAP analysis in- 
volves the assumption of a fiat universe, along with the strong 
HST prior on h and the age constraint to > 11 Gyr. For BBN 
the adopted primordial abundances are: yo = 10 5 (D/H) = 
2.6±0.4, Y(OSW) = 0.238±0.005, and Y(IT)= 0.244±0.002. 





Nv (2a bound) 


r/10 (2a range) 


WMAP 


8.3 


5.64 - 7.30 


yu + Y(OSW) 


3.3 


5.04 - 7.18 


yu + Y(IT) 


3.1 


4.89 - 6.56 


WMAP + y D + Y(OSW) 


3.3 


5.66 - 6.80 


WMAP + y u + Y(IT) 


3.1 


5.54 - 6.60 



TABLE II: The same as Table I, except that the constraint 
N v > 3 is imposed. 



Figs. 9 and 10. BBN and the primordial D abundance 
combine to provide a quite accurate determination of the 
baryon density (0.020 < w B < 0.025 at 2a). The cur- 
rently large uncertainty in the primordial abundance of 
4 He is responsible for the larger allowed range of N v . 
While the best fit value for N v is < 3 for either the OSW 
or the IT 4 He abundances, each are consistent with N v = 




FIG. 9: The marginalized likelihood distributions for 
ANv from the joint WMAP and BBN analysis for two choices 
of the primordial abundance of 4 He (solid: OSW, dashed: IT). 



9 



0.01 0.02 0.03 



/ 1 V 
1 / \l 

' ll 
1 I 

~~ ' y 

' 1^ 


11 

77 io ("b) : 


— r 

• 1 
ll 

fl 

ll 

ll 

|l 
ll 

V 

~ 2a 




l \ 



23456789 10 



FIG. 10: The marginalized likelihood distributions for 7710 
from the joint WMAP and BBN analysis for two choices of 
the primordial abundance of 4 He (solid: OSW, dashed: IT). 



3 at ~ 2a. These results are in excellent agreement with 
those of Hannestad [22], with which they are most di- 
rectly related. 

Clearly, BBN constrains N v much more stringently 
than WMAP, while the measurement of 7710 by WMAP 
is at a precision superior (by a factor ~ 2) to that from 
BBN. In this sense, the CBR and BBN are quite com- 
plimentary. Indeed, while the constraint on N v barely 



changes with the inclusion of the WMAP data in a joint 
analysis with BBN, it is sensitive to the adopted 4 He 
abundance; see Fig. 9. On the other hand, the joint con- 
straint on 77 is extremely insensitive to the choice of the 
4 He abundance, being dominated by the WMAP data 
(and the primordial D abundance); see Fig. 10. How- 
ever, BBN and WMAP do provide a very important 
consistency check of the standard model of cosmology 
at widely separated epochs, using significantly different 
physics. The excellent agreement between N v and 7710 
when the universe was 20 minutes and 380,000 years old 
is a major triumph for the (new) standard model of cos- 
mology. 
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